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Abstract 



We have developed a nonlinear differential effective dipole approximation 
(NDEDA), in an attempt to investigate the effective linear and third-order 
nonlinear susceptibility of composite media in which graded spherical inclu- 
sions with weak nonlinearity are randomly embedded in a linear host medium. 
Alternatively, based on a first-principles approach, we derived exactly the lin- 
ear local field inside the graded particles having power-law dielectric grada- 
tion profiles. As a result, we obtain also the effective linear dielectric con- 
stant and third-order nonlinear susceptibility. Excellent agreement between 
the two methods is numerically demonstrated. As an application, we apply 
the NDEDA to investigate the surface plasma resonant effect on the optical 
absorption, optical nonlinearity enhancement, and figure of merit of metal- 
dielectric composites. It is found that the presence of gradation in metal 
particles yields a broad resonant band in the optical region, and further en- 
hances the figure of merit. 

PACS number(s): 77.22.Ej, 42.65.-k, 42.79.Ry, 77.84.Lf 



1 



I. INTRODUCTION 



Graded materials, whose material properties can vary continuously in space, are abun- 
dant in nature. These materials have attracted much interest as one of the advanced in- 
homogeneous composite materials in various engineering applications [1]. With the advent 
of fabrication techniques, these materials can be well produced to tailor their properties 
for specific needs via the design of the material and microstructure gradients. Such a de- 
sign makes graded materials quite different in physical properties from the homogeneous 
materials and other conventional composite materials. Moreover, the composite media con- 
sisting of graded inclusions can be more useful and interesting than those of homogeneous 
inclusions. Although various theories have been established to investigate the optical and 
dielectric properties of the composite media of homogeneous inclusions [2,3], they fail to 
deal with the inhomogeneous composites of graded inclusions. Recently, a first-principles 
approach [4,5] and a differential effective dipole approximation [6,7] have been presented in 
order to investigate the dielectric response of graded materials. 

The problem becomes more complicated by the presence of nonlinearity in realistic com- 
posites. Besides inhomogeneity, such nonlinearity plays also an important role in the effec- 
tive material properties of composite media [8-11]. It is thus necessary to establish a new 
theory to study the effective nonlinear properties of graded composite media. In fact, the 
introduction of dielectric gradation profiles in nonlinear composites is able to provide an 
alternative way to control the local field fluctuation, and hence let us obtain the desired 
effective nonlinear response. 

In fact, the previous one-shell model [12] and multi-shell model [13], which were used to 
study the effective nonlinear optical property, can be seen as an initial model of graded inclu- 
sions. In this paper, we will put forth a nonlinear differential effective dipole approximation 
(NDEDA) to investigate the effective linear and nonlinear dielectric properties of composite 
media containing nonlinear graded spherical particles (inclusions). For such particles, the 
linear and nonlinear physical properties will continuously vary along their radius. 
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The paper is organized as follows. In section II, we describe the model and define briefiy 
the effective linear dielectric constant and third-order nonlinear susceptibility. In section III, 
the NDEDA is presented to investigate the effective linear dielectric constant and third-order 
nonlinear susceptibility of nonlinear graded composite media in the dilute limit. In section 
IV, based on a first-principles approach, we derive the exact solutions for composite media 
having power-law gradation profiles inside the inclusions, which is followed by the numerical 
results in section V. Finally, some conclusion and discussion is shown in Section VI. 

II. MODEL AND DEFINITION OF EFFECTIVE LINEAR AND NONLINEAR 

RESPONSES 

Let us consider a nonlinear composite system, in which identical graded spherical inclu- 
sions with radius a, are randomly embedded in a linear host medium of dielectric constant 
62- The local constitutive relation between the displacement (D) and the electric field (E) 
inside the graded particle is given by 



where e(r) and x(r) are, respectively, the linear dielectric constant and third-order nonlinear 
susceptibility. Note both e(r) and x(r) are radial functions. Here we assume that the 
weak nonlinearity condition is satisfied [14] . In other words, the contribution of the second 
(nonlinear) part [x.,(r)|Ep] in the right-hand side of Eq. (1) is much less than that of the 
first (linear) part e(r). We restrict further our discussion to the quasi-static approximation, 
under which the whole composite medium can be regarded as an effective homogeneous one 
with effective linear dielectric constant Cg and effective third-order nonlinear susceptibility 
Xe- To show the definitions of Ce and Xe, we have [14] 



where (• ■ ■) represents the spatial average, and Eq = EqBz is the external applied field along 



D = e(r)E + x(r)|E|2E, 



(1) 



(D) = eeEo + Xe|EorEo, 



(2) 



z axis. 
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The effective linear dielectric constant is given by 

eeEo = :^ 1^ eiEnn,idl/ = /(e(r)Eu„,i) + (1 - /)e2(Eu„,2), (3) 

where / is the volume fraction of the graded particles and the subscript stands for the linear 
local field [i.e., obtained for the same system but with x('^) = 0]. 

In view of the existence of nonlinearity inside the graded particles, Xe can then be written 
as [14,15] 

XeK = \ /^Xi|E|L,iEL,idl^ = \ x(r)|E|L,,EL,idy = /(x(r)|E|L,iEL,i). (4) 

In the next section, we will develop a NDEDA (nonlinear differential effective dipole 
approximation), in an attempt to derive the equivalent linear dielectric constant e(a) and 
third-order nonlinear susceptibility x(a) of the nonlinear graded inclusions. Then, the ef- 
fective linear dielectric constant and third-order nonlinear susceptibility of the composite 
media of nonlinear graded inclusions will be derived accordingly in the dilute limit. 

III. NONLINEAR DIFFERENTIAL EFFECTIVE DIPOLE APPROXIMATION 

To estabhsh the NDEDA, we first mimic the gradation profile by a multi-shell constric- 
tion. That is, we build up the dielectric profile by adding shells gradually [6] . We start with 
an infinitesimal spherical core with linear dielectric constant e(0) and third-order nonlinear 
susceptibility x(0), and keep on adding spherical shells with linear dielectric constant e(r) 
and third-order nonlinear susceptibility xi''^) at radius r, until r = a is reached. At radius r, 
the inhomogeneous spherical particle with space-dependent dielectric gradation profiles e(r) 
and x(r) can be replaced by a homogenous sphere with the equivalent dielectric properties 
e(r) and xi^)- Here the homogeneous sphere should induce the same dipole moment as the 
original inhomogeneous sphere. 

Next, we add to the sphere a spherical shell of infinitesimal thickness dr, with dielec- 
tric constant e(r) and nonlinear susceptibility x('")- this sense, the coated inclusions is 
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composed of a spherical core with radius r, hnear dielectric constant e(r) and nonlinear 
susceptibility ^ shell with outermost radius r + dr, linear dielectric constant e(r) 

and nonlinear susceptibility x(r). Since these coated inclusions are randomly embedded in a 
linear host medium, under the quasi-static approximation, we can readily obtain the linear 
electric potentials in the core, shell and host medium by solving the Laplace equation [16] 

0c = -EqAR cos 6, R <r, 

(f)s^ -Eoi BR - ^ \ cose, r<R<r + dr, (5) 

/ D(r + drf\ 

(Ph = -EoiR j ^Q^^' R>r + dr, 

where 

962e(r) ^ 362[e(r) + 26(r)] ^ _ 3e2[e(r) - e(r)] 



D 



g ' Q Q 

[e(r) - e2][e(r) + 2e(r)] + + 2e(r)][e(r) - e(r)] 



with interfacial parameter A = [r/(r + dr)]^, and 

g = [e(r) + 262] [e(r) + 2e(r)] + 2A[e(r) - e2][e(r) - e(r)]. 

The effective (overall) linear dielectric constant of the system is determined by the dilute- 
limit expression [17] 

ee = €2 + 3pe2-D, (6) 

where p is the volume fraction of graded particles with radius r. The equivalent dielectric 
constant e{r+dr) for the graded particles with radius r+dr can be obtained self-consistently 
by the vanishing of the dipole factor D by replacing €2 with e{r+dr). Taking the limit dr — > 
and keeping to the first order in dr, we obtain 

e(r) - e(r) 



e(r -|- dr) = e(r) -|- 3e(r)A 



e(r)(l-A)+e(r)(2 + A) 
eir) — e(r 



3 + 



e(r) 



dr. (7) 



Thus, we have the differential equation for the equivalent dielectric constant e(r) as [6] 

de(r) _ [e(r) - e(r)] ■ [e(r) + 2e(r)] 



(8) 



dr re(r) 

Note that Eq. (8) is just the Tartar formula, derived for assemblages of spheres with varying 
radial and tangential conductivity [18]. 

Next, we speculate on how to derive the equivalent nonlinear susceptibility x(r). After 
applying Eq. (4) to the coated particles with radius r + dr, we have 



x(r + dr) 



(|EpE^)R<r+dr 



(|E 


|2-p2\ 
1 -CJ jR<r 




Eq 


2E§ 



^ + (1 - A) 



r<R<r+Ar 

WEg ■ 



(9) 



As dr ^ 0, the left-hand side of the above equation admits 



x(r + dr) 



(|E| 


^ E^)ij<r+dr 




|Eo 


mi 



e(r + dr) + 2e2 



3eg 



x{r)\K\^K^-drxir)\K\^K^ 
dr. 



V 

^e{r + dr) + 262 J 
3de(r)/dr ^ f de(r)/dr 



2e2 + e(r) \2e2 + e{r) 



dr 



(10) 



with K = (862)/ [e(r) + 2e2]. The first part of the right-hand side of Eq. (9) is written as 



|E„|2E5 



dr 



l + (6y + 2y*-3) — 



(11) 



where 



y 



[6(r)-62][e-(r)-e(r)] 



e(r)[e(r) + 2e2] " 
The second part of the right-hand side of Eq. (9) has the form [17] 

^ ^ lEoPEg 5r ' ' 



x(5 + I8x^ + 18|a;|2 + 4a;^ + 12a;|a;|^ + 24|a;|V), (12) 



where 



''':'-'''•> and ,= ^^|fM + 2^M] 



e(r) + 2e(r) " " e(r)[e(r) + 2e2j 
Substituting Eqs. (10), (11) and (12) into Eq. (9), we have a differential equation for the 

equivalent nonlinear susceptibility x(r), namely. 
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dx(r) 
dr 



3de(r)/dr ^ f de(r)/dr 



2e2 + e(r) V2e2 + e(r), 



+ X(r) • 



6y + 2y* 



3 , 3x(r) 
5r 



e(r) + 2e(r) 



3e(r) 



^e(r) + 2e(r) 
. 3e(r) 



(5 + 18x^ + 18|xp + 4x'^ + 12x\x\^ + 24|a;| V) 



(13) 



So far, the equivalent e(r) and x(^) of graded spherical particles of radius r can be cal- 
culated, at least numerically, by solving the differential equations Eqs. (8) and (13), as long 
as e(r) (dielectric-constant gradation profile) and x('") (nonlinear-susceptibility gradation 
profile) are given. Here wc would like to mention that, even though xi^^) is independent 
of r, the equivalent x(?^) should still be dependent on r because of e(r) as a function of r. 
Moreover, for both e(r) = ei and x('") — Xi (i-6-, they are both constant and independent 
of r), Eqs. (8) and (13) will naturally reduce to the solutions e(r) = ei and x{r) = Xi- 

To obtain e(r = a) and x{i^ — o), we integrate Eqs. (8) and (13) numerically at given 
initial conditions e(r — > 0) and xir — > 0). Once e(r = a) and xir — a) are calculated, we 
can take one step forward to work out the effective linear and nonlinear responses eg and Xe 
of the whole composite in the dilute limit, i.e. [14], 

e(r = a) - 62 



ee^e2 + 362 f 



e(r ^ a) + 2e2 ' 



(14) 



and 



Xe = fxir = a) 



3e2 

e(r ^ a) + 2e2 



e r 



3e2 

= a) + 262 



(15) 



IV. EXACT SOLUTION FOR POWER-LAW GRADATION PROFILES 

Based on the first-principles approach, we have found that, for a power-law dielectric 
gradation profile, i.e., e(r) = A(r/a)", the potential in the graded inclusions and the host 
medium can be exactly given by [4] 

0i(^) — —CiEgr'^ cos9, r < a, 

(f)h{r) ^ -Eorcos9+^Eocos9, r > a, (16) 
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where the coefficients and ^2 have the form 



and s is given by 



SA + 262 



1 



sA-e2 3 
and to — a , 

^ sA + 262 



'=2 



V9 + 2n + - (1 + n) 



The local electric field inside the graded inclusions can be derived from the potential 
E = -V(/), 

Ej = ^iEQr^~^{s cos der — sin^e^) = {i-E'or^~^{(s — 1) cos ^ sin ^ cos 

+ (s - I)cos6'sin6'sin0ej/+ [{s - 1) cos^ 6' + Ije^}, (17) 

where e^, e^, e^, and are unix vectors. In the dilute limit, from Eq. (3), we can obtain 
the effective hnear dielectric constant as follows 



£2 + / [A{r/aY - £2] e, • E,dy 



62 + 362/- 



2 + s 



£2 



sA + 2€2\2 + n + s 2 + s 
On the other hand, the substitution of Eq. (17) into Eq. (4) yields 

1 



(18) 



Xe^TT f x(r)|6rC?(s^cos2^ + sin2^)V*-2sinMrdM0 

Jo 



(19) 



For example, for a linear profile of x(r), i.e., x(r) = /ci + /c2 ■ r/a, Eq. (19) leads to 



'(8 + 4.^ + 3/) 



/ 


3e2 


2 


f 362 \ 


20 


sA + 262 




KsA + 262) 



As -I 



(20) 



In addition, for a power-law profile of x(r), namely, x(^) = ki{r/a)'''^, Eq. (19) produces 



Xe 



/ 



3e2 



sA + 262 



362 / 8 + 4g^ + 3g^ ' 

sA + 262/ H ^2 - 1 + 4s 



(21) 
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V. NUMERICAL RESULTS 



We are now in a position to evaluate the NDEDA. For the comparison between the 
first-principles approach and the NDEDA, we first perform numerical calculations for the 
case where the dielectric constant exhibits power-law gradation profiles e(r) = A{r/a)^, 
while the third-order nonlinear susceptibility shows two model gradation profiles: (a) linear 
profile x(r) = ki + k2 ■ r/a, and (b) power-law profile x(^) = ki{r/a)^^. Without loss of 
generality, we take £2 = 1 and a = 1 for numerical calculations. The fourth-order Runge- 
Kutta algorithm is adopted to integrate the differential equations [Eqs. (8) and (13)] with 
step size 0.01. Meanwhile, the initial core radius is set to be 0.001. It was verified that this 
step size guarantees accurate numerics. 

In Fig. 1, the effective hnear dielectric constant (ee) is plotted as a function of A for 
various indices n. It is shown that eg exhibits a monotonic increase for increasing A (and 
decreasing n). This can be understood by using the equivalent dielectric constant e(r = a) 
which increases as A increases {n decreases). Moreover, the excellent agreement between 
the NDEDA [Eq. (8)] and the first-principles approach [Eq. (18)] is shown as well. 

Next, the effective third-order nonlinear susceptibility (xe) is plotted as a function of A 
for the linear gradation profile x(r) — ki + k2 ■ r/a (Fig. 2), and for the power-law profile 
x(r) = ki{r/aY'^ (Fig. 3). We find that the effective nonhnear susceptibility decreases for 
increasing A. The reason is that, as mentioned above, for larger A, the graded inclusions 
possess larger equivalent dielectric constant, and the local field inside the nonlinear inclusions 
will become more weak, which results in a weaker effective nonlinear susceptibility (xe)- In 
addition, increasing n leads generally to increasing Xei and such a trend is clearly observed 
at large A. Again, we obtain the excellent agreement between the first-principles approach 
[Eqs. (20) and (21)] and the NDEDA [Eqs. (8) and (13)]. 

In what follows, we investigate the surface plasma resonance effect on the metal-dielectric 
composite. We adopt the Drude-hke dielectric constant for graded metal particles, namely. 
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where uip{r) and 7(r) are the radius-dependent plasma frequency and damping coefficient, 
respectively. For the sake of simplicity, set x(?^) = Xi to be independent of r, in an attempt to 
emphasize the enhancement of the effective optical nonlinearity, and 62 = 1.77 (the dielectric 
constant of water). We assume further ujp{r) to be 

T 

u!p{r) = (jjp{l - ku, ■ -) , r < a. (23) 

(X 

This form is quite physical for k^^ > 0, since the center of grains can be better metallic so 
that ujp{r) is larger, while the boundary of the grain may be poorer metallic so that ujp{r) 
is much smaller. Such the variation can also appear because of the temperature effect [19]. 
For small particles, we have the radius-dependent 7(r) as [20] 

k 

7{r) = 7(00) + r<a, (24) 

r/a 

where 7(cxo) stands for the damping coefficient in the bulk material. Here k^ is a constant 
which is related to the Fermi velocity vp- In this case, the exact solution being predicted by 
a first-principles approach is absent. Fortunately, we can resort to the NDEDA instead. 

In Fig. 3, we plot the optical absorption [~ Im(ee)], the modulus of the effective third- 
order optical nonlinearity enhancement (|Xe|/xi) and the figure of merit (|xe|/Im(ee)) versus 
the incident angular frequency cu. For various variance slopes ki_^, large surface plasma reso- 
nance peaks can appear around co ~ OAup, as expected. However, when the radius depen- 
dence of a;p(r) is taken into account (i.e., ki_^ ^ 0), besides a sharp peak, a broad continuous 
resonant bands in the high frequency region is clearly observed. The broad spectra result 
from the effect of radius dependence of the plasma frequency. In Ref. [12], we found that, 
when the shell model is considered, there should be a broad continuous spectrum around the 
large pole in the spectral density function. In fact, the graded particles under consideration 
can be regarded as a certain limit of multi shells, which hence should lead to broader spectra 
in Im(ee), |Xe|/Xi as well as |xe|/Ini(ee)- In addition, we note that, with increasing ki_j, both 
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the surface plasma frequency and the center of resonant bands are red-shifted. In particular, 
the resonant bands can become more broad due to strong inhomogeneity of the particles. 
From the figure, we conclude that, although the third-order optical nonlinearity is always 
accompanied with the optical absorption, the figure of merit in the high frequency region is 
still attractive due to the presence of weak optical absorption. Thus, we believe that graded 
particles may have potential applications in obtaining the optimal figure of merit, and make 
the composite media more realistic for practical applications. 

Finally, we focus on the effect of 7(r) on the nonlinear optical property in Fig. 5. As 
evident from the results, the variation of plays an important role in the magnitude of the 
effective optical properties, particularly at the surface plasma resonance frequency. 

VI. CONCLUSION AND DISCUSSION 

Here a few comments are in order. In this work, we have developed an NDEDA (nonlin- 
ear differential effective dipole approximation) to calculate the effective linear and nonlinear 
dielectric responses of composite media containing nonlinear graded inclusions. The re- 
sults obtained from the NDEDA are compared with the exact solutions derived from a 
first-principles approach for the power-law dielectric gradation profiles, and the excellent 
agreement between them has been shown. We should remark that the exact solutions are 
also obtainable for the linear dielectric gradation profiles with small slopes (the derivation 
not shown here). In this case, the excellent agreement between the two methods can be 
shown as well since the NDEDA is vahd indeed for arbitrary gradation profiles. In general, 
the exact solution is quite few in realistic composite research, and thus our NDEDA can be 
used as a benchmark. 

It is instructive to develop the first-principles approach to nonlinear graded composites. 
The perturbation approach [21] in weakly nonlinear composites as well as the variational 
approach [22] in strongly nonlinear composites are just suitable for this problem. 

The NDEDA is strictly valid in the dilute limit. To achieve the strong optical nonlinearity 
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enhancement, we need possibly nonlinear inclusions with high volume fractions. In this 
connection, the effect of the volume fraction is expected to cause a further broadening of 
the resonant peak, and possibly, a desired separation of the optical absorption peak from 
the nonlinearity enhancement due to mutual interactions [23]. Therefore, it is of particular 
interest to generalize the NDEDA for treating the case of high volume fractions. 

To one's interest, the NDEDA can be applied to biological cells. Since the interior 
of biological cells is often inhomogeneous and nonlinear in nature, they can be treated as 
particles having dielectric gradation profiles [7]. Moreover, in the case of cell membranes 
containing mobile charges introduced by the adsorbed hydrophobic ions, the local dielectric 
anisotropy occurs naturally, and should be expected to play a role [24]. Thus, it is also 
interesting to see what happens to the NDEDA as one takes into account the local dielectric 
anisotropy. The resultant anisotropic NDEDA will help to investigate the AC electrokinetic 
phenomena of biological cells [25]. Work along this direction is in progress. 

To sum up, we put forth an NDEDA and a first-principles approach for investigating the 
optical responses of nonlinear graded spherical particles. The excellent agreement between 
the two methods has been shown. As an apphcation, we applied the NDEDA to discuss 
the surface plasma resonance effect on the effective linear and nonlinear optical properties 
like the optical absorption, the optical nonlinearity enhancement, and the figure of merit. 
It is found that the dielectric gradation profile can be used to control the surface plasma 
resonance and achieve the large figure of merit in the high-frequency region, where the 
optical absorption is quite small. 
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FIGURES 

FIG. 1. The effective linear dielectric constant eg versus A for the power- law dielectric gradation 

profile e(r) = A{r/a)" in the dilute limit / = 0.05. Lines: numerical results from the NDEDA 
[Eq. (8)]; Symbols: exact results [Eq. (18)]. 

FIG. 2. The effective third-order nonlinear susceptibility Xe versus A for power-law dielec- 
tric-constant gradation profile e(r) = A(r/a)" and linear nonlinear-susceptibility gradation profile 
x(r) = ki + k2 ■ r/a with (a) ki = 1 and ^2 = 1, and (b) ki = 2 and ^2 = 3. Lines: numerical 
results from the NDEDA [Eqs. (8) and (13)]; Symbols: exact results [Eq. (19)]. 

FIG. 3. Same as Fig.2, but for power-law nonlinear-susceptibility gradation profile 
x(r) = ki{r/a)^'^. 

FIG. 4. (a) The linear optical absorption Im(ee), (b) the enhancement of the third-order op- 
tical nonlinearity |xe|/xi) and (c) the figure of merit = |xe|/Im(ee) versus the incident angu- 
lar frequency u/tOp for dielectric-constant gradation profile e(r) = 1 — LOp{r)/[uj{LO + ilif))] with 
ujp{r) = u!p{l — k^ ■ r/a) and j{r) = 0.01a;p. Parameters: e2 = 1.77 and / = 0.05. 

FIG. 5. Same as Fig.4, but with oJp{r) = ujp and j{r) = 7(00) -|- k^/{r/a) for ^{oo) = O.OloOp. 
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Fig. 2. /Gao, Huang and Yu 
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Fig.3. /Gao, Huang and Yu 
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